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ABSTRACT 

Mesh adaptation for finite element approximation is a proce- 
dure used in numerous applications. The use of thin and long 
anisotropic triangles improves the efficiency of the procedure. 

When piecewise linear finite elements are used, the aspect 
ratio for mesh adaptation is generally dictated by the absolute 
value of the (estimated) hessian matrix of the approximated 
function. We give in this paper the corresponding aspect ratio 
for piecewise quadratic finite elements. 

Keywords — Anisotropic finite elements, Adaptive meshes, 
Interpolation, Nonlinear approximation. 

1. INTRODUCTION 

Consider a bounded polygonal domain C R^, a sufficiently 
smooth function / : — > R, and an integer m > 2. We 
introduce the problem of optimal mesh adaptation 

min{#(r); Ts.t. || V(/ - I^^" V)IU^(o) < e}, (D 

where T stands for an arbitrary triangulation of il, and #(T) 
for its cardinality. Here I™^^ denotes the Lagrange interpola- 
tion operator onto finite elements of degree m — 1 on T. 

In practical applications, the problem ([T]) is generally in- 
tractable for at least three reasons. 1: The function / may 
have complicated local features, difficult to analyze. We thus 
first make a local analysis based on Taylor developments. 2: 
The collection of triangular meshes of 51 is a combinatorial set 
and problems such as ([T]) are typically NP-complete (after dis- 
cretization). We avoid this problem by first considering the case 
of a single triangle. 3: Currently available anisotropic mesh 
generation algorithms only give control on the aspect ratio and 
orientation of the generated triangles, but not on their other fea- 
tures. We thus only optimize this aspect ratio. 

2. AN OPTIMIZATION PROBLEM 

We denote by Pm-i the space of bivariate polynomials of de- 
gree < m — 1, and by H„j the space of homogeneous polynomi- 
als of degree m. If / G C™(f2), if z G is fixed and if /i G R^ 
is small, then locally 



for some fiz S Pm-i and tt^ € H„j. If T is a sufficiently small 
triangle, we thus have at least heuristically on T 

V(/-irV)^V(7r,-I^-i7r,), (3) 

since the Lagrange interpolation operator I™~^ on the triangle 
T reproduces the elements of P,n-i. 

For any triangle T and any / e H^{T) n C"(T), we define 
the averaged interpolation error eT{f)m as follows 

eriffm :-^^|V(/-I??-V)P- 

The local counterpart of ([TJ is the problem of the optimal trian- 
gle : find for all tt G H,„ 

sup{|r|; Ts.t. ctW™ < 1}. (4) 

Indeed the cardinality of a triangulation is inversely propor- 
tional to the area of its elements. This approach is developed 
in Chapter 2 of and leads to asymptotically optimal error 
estimates of ([T]l as e — > (or more precisely estimates of e 
as #(T) — > oo, which is equivalent). Unfortunately these 
estimates are not completely realistic for applications because 
currently available numerical anisotropic mesh generators only 
control the aspect ratio and orientation of the generated trian- 
gles. 

For each triangle T, of vertices vi, V2 and v^, we denote by 
zt := (wi + U2 + W3)/3 its barycenter We denote by S*^ the 
collection of 2 x 2 symmetric positive definite matrices, and we 
define a matrix T-Lt G 5*^ by the equality 

2 

l<i<3 

If A is an invertible 2x2 matrix and if T' is mapped onto T by 
the linear map z M> Az, then one easily checks that 

Ut' = A^HtA. (5) 

By construction the triangle Tcq of vertices 
(cos(2fc7r/3), sin(2fc7r/3))o<fc<2 satisfies Ht^^ = Id. 
Combining these two properties. Proposition 5.1.3 in 111 
establishes that for any triangle T 



f{z + h)^ fi,{h)-^n,{h)-^o{\hr), 



(2) 



|r|VdetHT = |Teq|, 





Fig. 1. Lagrange interpolation points for Pi and P2 finite ele- 
ments (left), triangle T and associated ellipse St (right). 

and that there exists a rotation U (depending on T) such that 

z ^ UH^z - zt) (6) 

maps T onto T^q (the power a of a symmetric positive definite 
matrix is obtained by elevating the eigenvalues to the power a in 
a diagonalization). Furthermore the ellipse of minimal volume 
containing T is £t {z; (2 — zt)^'Ht{z — zt)}, see Fig 1. 
The matrix Ht thus encodes the area, the aspect ratio and the 
orientation of T. 

For each M e S2 and each tt e we define 

eAf(7r)m := sup{eT(7r)m; T s.t. Ht = A/}. 

We finally introduce for each tt e H„ the problem of the opti- 
mal aspect ratio for Pm-i interpolation 



inf{detM; M e s.t. eji/(7r)™ < 1}. 



3. MAIN RESULT 



(7) 



Our main result is the solution of the optimization problem (j7]i 
in the case of piecewise linear and piecewise quadratic finite el- 
ements. The piecewise quadratic case is entirely new and gives 
a well founded answer to a long standing question: which as- 
pect ratio, depending on the third derivatives of the approxi- 
mated function, should be used in finite element software that 
combine anisotropy and P2 elements ? 

We first introduce some notation. We equip the vector space 
H„i with the norm 

Ikll sup |7r(M)|. 

|u|<l 

For each tt e H2, tt — ax^ + 2hxy + cy^, we define 

a b 
b c 

The absolute value of a symmetric matrix (resp. the square root 
of a non negative symmetric matrix) is obtained by taking the 
absolute value (resp. square root) of the eigenvalues in a diago- 
nalization. For each tt £ H2 we set 

M2{tt) :=||7r|| |[7r]|. 

For each tt e H3, tt = ax'^ + ibx^y + 3cxy^ + dy^, we set 

^ 1 
' disc TT , 

Id, 



MsiTT) := J[d,TT]^ + [dyTT]^ + 



where disc tt := 4(ac — b'^)(bd — c^) — {ad — bc)"^ and := 
max{A, 0}. 

Theorem. For m E {2,3} the map tt e — > A^,„(7r) is 
a near-minimizer of the problem Q in the following sense. If 
TT is non-univariate then Mrn{TT) is non-degenerate. Further- 
more there exists a constant C, independent of tt, such that 

eAl„(7r)(7r)m < C and 

detA^,„(7r) < Cinf{det M; M e 5^ s.t. eA/(7r),„ < 1}. 
Proof: The integer m e {2, 3} is fixed, and we denote for each 

TT EMm 

IIVttII sup |V7r(u)|. 

|«|<i 

For each 2x2 matrix A we denote hy tt o A the element of 
defined by (tt o A){u) := tt{A(u)), u e R^. We recall that 

V(7r o A){u) = A^Vtx{A{u)), u e R^^ 

which implies for any rotation U 

IIVttII = ||V(7ro C/)||. (8) 

The main difficulty of this proof is to show that there exists a 
cons 
has 



constant C > 1 such that for aU tt € and all M e 5*2 one 



C-^eM{TTU < ||M||5||V(^oAf-5)|| <CeM(7r), 



(9) 



Assume that this point is established. Proposition 6.5.4 in HI, 
states that the map tt !—> M„i{tt) is a near-minimizer for the 
optimization problem 

inf{detM; M G s.t. ||M|| 5 || V(7r o Af-5 )||}^ 

in the same sense as in the statement of this theorem. Com- 
bining this result with the equivalence (|9]), and using the homo- 
geneity of TT, we immediately conclude the proof of this theo- 
rem. 

We thus turn to the proof of (j9]l. Our first observation is that 
there exists a constant Cq such that for all tt € 



1 



cql 



iv(^-i":V)p<Coiiv^ii, 



(10) 

indeed the left and right hand side are norms on Hm. 

Consider a symmetric matrix M e S2 and a triangle T such 
that Ht = M. According to (|6]l there exists a rotation U such 
that the image of T by the map z 1— > UAI^^z — zt) is the 
triangle Tcq. Injecting this change of variables in ([TO]) we obtain 




vr)P<Co||V(^o(M-5[/-i))||. 



Observing that ||Az)|| > \\A ^|| ^|u| for any invertible 2x2 
matrix A and vector 1; G R^, and recalling M, we obtain 

||M||-5eTW™ <Co||V(^oM-5)||. 



Taking the supremum of the left hand side among all triangles 
T such that Ht = M we establish the left part of (j9]l, provided 
thatC > Co. 

We now remark that there exists a constant Ci such that for 

all TT e Hm 



1) 



Indeed assume that the right hand side vanishes. Then fi := 
TT — I™^^ TT is a polynomial of degree m depending only on the 
variable y, and which vanishes on the Lagrange interpolation 
points of Teq, see Figl. Hence ii vanishes for y — ±\/3/2 and 
y = if m = 2 (resp. y = ±V3/2, y = ±\/3/4 and ?; = if 
m — 3). Therefore /i = which implies that n — 0. Both sides 



of ( 1 1 1 are thus equivalent norms on the vector space . 
We consider a diagonalization of a symmetric matrix Ma, 

M e Si, 

' a 



M5 = U'^DU, D = 



/3 



where U is & rotation and a = 2 . Consider the triangle T 
which is mapped onto Tcq by the change of coordinates 

z ^ UM^z = DUz, 

and thus satisfies Ht = M according to Q. Injecting this 
change of variables into ( 1 1 1 we obtain 



||V(7ro(Af-5[/-i))|| < Ci 




where v := U ^Cx, e.^ :— (1, 0), and where we used for the dx 
derivative that J7^^_D^^ = a~^v. Recalling that a — ||M||2, 
\v\ = 1, and using ^ we obtain 



|M|l5||V(7roM"5)|| < C 



= CieT(7r)„i < CieM(7r)m. 



This concludes the proof of (j9| with C :~ maxjCo, Ci}, hence 
the proof of this theorem. o 

4. APPLICATIONS AND CONCLUSION 

Consider a function / for which one desires to solve, at least 
heuristic ally, the optimization problem ([T]l. Assume that some 
estimate of 



k+l—m 



d^xd'-y kl II 



is known at each point z ^ il, and define a riemannian metric 
on as follows 



H{z) := A(det7W™(7r^)) -^-^Mmiiiz), 



(12) 
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Fig. 2. Interpolation of ( 13 1 with Pi elements on a uniform. 



isotropic or anisotropic mesh of cardinality 500. 

where A > is a constant (this expression needs to be slightly 
modified if tt^ vanishes or is univariate for some values of 2, 
in order to ensure that H G (7*^(11, S"^)). Some mesh gen- 
erators such as |2| can, at least heuristically, and provided H 
has sufficient regularity, produce a mesh T of such that 
C-^H{z) < Ht < CH{z) for each T e T each z e T, 
where C is a constant not too large. In other words the aspect 
ratio of the elements of T is dictated by the metric H. Some 
rigorous results in this direction can be found in Chapter 5 of 

m. 

In the expression ( 12 1 the matrix M-miT^z) ensures that the 
elements of T have the optimal aspect ratio, while the scalar fac- 
tor (det A^m(7rz))~ 2^ guarantees that the interpolation error is 
equidistributed among the elements of T (a general principle in 
adaptive approximation). 

We conducted some numerical experiments using ||2l and for 
the synthetic function 



f{x,y) := tanh(10(sin(5y) - 2x)) 



x'^y - 



(13) 



on the domain fl (—1,1)-^. They illustrate the improvement 
offered by anisotropic mesh adaptation, both in the case of Pi 
and P2 elements, for a triangulation of cardinality 500. 



#(r) = 500 
|v(/-i^/)||l2 
|v(/-i^/)||l2 



Uniform 
110 
79 



Isotropic 
51 
14 



anisotropic 
11 
0.88 



Our next objective is to combine our analysis with an adaptive 
anisotropic mesh refinement procedure, for a partial differential 
equation solved with P2 finite elements. The optimization prob- 
lem ([T]) is particularly relevant in the case of elliptic equations. 
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